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We present a rigorous study of vibrational relaxation in p-Hb + p-Eb collisions at cold and ul- 
tracold temperatures and identify an efficient mechanism of ro-vibrational energy transfer. If the 
colliding molecules are in different rotational and vibrational levels, the internal energy may be 
transferred between the molecules through an extremely state-selective process involving simulta- 
neous conservation of internal energy and total rotational angular momentum. The same transition 
in collisions of distinguishable molecules corresponds to the rotational energy transfer from one 
vibrational state of the colliding molecules to another. 



The creation of ultracold molecules has opened up new 
avenues for research in physics and chemistry. Ultracold 
molecules are finding exciting applications in condensed- 
matter physics [H, quantum information science Q, 
physics of elementary particles [H and chemical reaction 
dynamics [ij. The experimental work on the production 
of dense ensembles of ultracold molecules is currently be- 
ing pursued by many research groups The progress 
of the experiments relies on detailed understanding of 
molecular collisions at ultracold temperatures [f| @ . Re- 
cent theoretical studies @, H, @j 03 have provided impor- 
tant insights into atom - molecule inelastic and reactive 
collisions at cold and ultracold temperatures. The colli- 
sion energy of ultracold molecules is usually much smaller 
than the energy spacing between ro-vibrational molecu- 
lar levels and the duration of an ultracold collision is very 
long, which leads to unusual propensity rules and mech- 
anisms of inelastic scattering. For example, vibrational 
relaxation in ultracold collisions of molecules with atoms 
was found to give rise to quasi-resonant energy trans- 
fer [TTI | resulting in extremely narrow distributions over 
final ro-vibrational states. Ultracold molecules in excited 
vibrational states are routinely produced in experiments 
on photoassociation and Feshbach resonance linking of 
ultracold atoms [12] • The lifetime of the molecules is lim- 
ited by collision-induced vibrational relaxation. Dense 
ensembles of ultracold molecules can also be created in 
experiments on buffer gas cooling [l| or Stark deceler- 
ation of molecular beams [l4j. Ground-state molecules 
may then be excited by black-body radiation fl5j and 
release energy through vibrationally inelastic collisions. 
Laser cooling of molecules, when realized, will result 
in significant populations of vibrationally excited states 
that may relax through collisions in a gas of sufficient 
density [161 ]. Evaporative cooling of molecular ensem- 
bles to ultracold temperatures is based on energy dissi- 
pation in molecule - molecule collisions. The spectrum 
of accessible energy levels is much denser for collisions of 
molecules with each other than for collisions of molecules 
with atoms and it is important to understand dominant 
energy transfer mechanisms in molecule - molecule inter- 
actions at ultracold temperatures. 



In this Letter, we present rigorous quantum calcula- 
tions of cross sections for cold and ultracold collisions 
of p-H2 molecules in different levels of vibrational and 
rotational excitation and identify a new mechanism for 
ro-vibrational energy transfer, not present in collisions of 
molecules with atoms. We find that the energy exchange 
process in collisions of indistinguishable molecules at cold 
and ultracold temperatures may be highly efficient and 
state selective when the internal energy and total rota- 
tional angular momentum of the colliding molecules are 
simultaneously conserved. By treating the H2 molecules 
as distinguishable particles it is shown that this mecha- 
nism is primarily driven by a rotational transfer process. 
We show that the efficiency of the near-resonant energy 
transfer is largely insensitive to the initial vibrational ex- 
citation of the molecules. 

The H2 + H2 system is the simplest neutral tetra- 
atomic system and it serves as a prototype for describing 
collisions between diatomic molecules. Using a model 
rigid rotor potential energy surface (PES) for the H2-H2 
system developed by Zarur and Rabitz [l?!], Forrey [ill ] 
investigated rotational transitions in P-H2 + P-H2 col- 
lisions at ultracold temperatures. Boothroyd, Martin, 
Keogh and Peterson (BMKP) have recently calcu- 
lated a full-dimensional PES for the H4 system. Pogreb- 
nya and Clary [20( reported rate coefficients for collisions 
of P-H2 and 0-H2 in the first excited vibrational state 
for collision energies between 1(T 3 eV (11.604 K) and 
1 eV (11604 K). They included all internal degrees of 
freedom and used the BMKP PES in a quantum calcula- 
tion based on the centrifugal decoupling approximation. 
Mate et al. [2l| reported an experimental study of rota- 
tionally inelastic collisions at temperatures between 2 K 
and 110 K and compared their data to theoretical results 
obtained using the PESs calculated by Diep and John- 
son (DJ) [H and Schaefer and Kohler (SK) These 
potentials describe the interaction between H2 molecules 
fixed at the equilibrium internuclear distance. Lee et 
al. [24j have recently presented a comparative analysis of 
cross sections for rotationally inelastic collisions between 
P-H2 at low and ultralow energies computed using the D J 
and BMKP PESs. Time-dependent approaches have also 
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been recently used to describe the dynamics between two 
P-H2 [HI, Hgj and between p-H 2 and o-H 2 [27| at thermal 
energies. 

The results presented here are obtained using the 
BMKP PES and a full-dimensional quantum dynam- 
ics calculation without any decoupling approximations. 
We employed a newly developed quantum scattering 
code [28| based on the formalism described by Arthurs 
and Dalgarno 29], Takayjinagi [30| . Green [31] • and 
Alexander and DePristo [321. We consider H2 molecules 
initially in different quantum states characterized by the 
vibrational quantum number v and the rotational an- 
gular momentum j. We refer to a combination of two 
ro-vibrational states of H 2 as combined molecular state 
(CMS). A CMS represents a unique quantum state of the 
diatom - diatom system before or after a collision. The 
CMS will be denoted as (vjv'j'). We used a "well-ordered 
states" notation [3(3] where v > v' or when v = v', 
j > j' . As the colliding molecules are indistinguishable, 
the quantum numbers vj or v'j' cannot be assigned to a 
particular molecule. To describe collisions of distinguish- 
able molecules, we will use the notation [vj;v'j'], where 
the quantum numbers vj characterize one molecule, and 
v'j' the other. The cross sections for state-resolved ro- 
vibrational transitions are obtained from the solution of 
coupled differential equations (3lj | using the log-derivative 
method of Johnson [33| and Manolopoulos |34| , and ap- 
propriate scattering boundary conditions [29|. 

Fig. [Upresents the cross sections for elastic and inelas- 
tic scattering of H 2 molecules for collision energies from 
1CP 6 K to 120 K. We consider collisions of H 2 molecules 
in the excited ro-vibrational states (v = l,j = 0) and 
(v = l,j = 2) with H 2 molecules in (v = 0, j = 0) and 
collisions of H 2 (u = l,j = 0) with H 2 (w = 0,j = 2). The 
cross sections for elastic collisions of H 2 molecules in the 
ground state (v = 0, j = 0) are also shown. The cal- 
culations are performed for the total angular momenta 
J = [0 — 10] and for both parities. While the elastic 
scattering cross sections for collisions of molecules in the 
vibrationally excited states are practically the same, the 
probability of inelastic scattering is extremely sensitive to 
the initial states of the colliding molecules. The inelastic 
relaxation of H 2 (v = 1 , j = 0) is six orders of magnitude 
more efficient in collisions with H 2 (i> = 0, j = 2) than in 
collisions with H 2 (t> = 0, j = 0). The zero-temperature 
limiting values of the relaxation rate constants for (1000), 
(1200) and (1002), are respectively, 9 x 10~ 18 , 4 x 10~ 14 
and 6 x 10~ 12 cm 3 s _1 . The inelastic cross section for 
collisions between H 2 (v = \,j = 2) and H 2 (t> =0,j = 0) 
is smaller than the inelastic cross section for collision be- 
tween H.2(v = l,j = 0) and H 2 (t> = 0,j = 2) at low and 
ultralow energies. It undergoes a rapid enhancement of 
two orders of magnitude at the collision energy 25.45 K, 
at which both inelastic cross sections become similar. 

To explain the results of Fig. [I] we have calculated the 
corresponding state-to-state cross sections for a collision 
energy of 10~ 6 K (see Fig. [2]). The final CMSs are pre- 
sented in the order of increasing energy. Figure [5] shows 
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FIG. 1: Elastic (upper panel) and inelastic (lower panel) cross 
sections for H2-H2 collisions for collision energies from 1(T 6 K 
to 120 K. The curves are labeled by {vjv'j') values. See text 
for details. 



that ro-vibrational relaxation in molecule - molecule col- 
lisions is determined by the propensity to conserve the 
internal energy and the total rotational angular momen- 
tum of the colliding molecules. The cross sections for vi- 
brational relaxation in collisions of H 2 (u = l,j = 0) with 
H 2 (u — 0,j — 0) are rather small because the minimiza- 
tion of the energy gap between the initial and final CMS 
would require a large change of the total rotational an- 
gular momentum of the colliding molecules. Conversely, 
the conservation of the total rotational angular momen- 
tum would introduce a large energy gap. Consequently, 
there is no single dominant CMS in the final state-to- 
state distribution shown in panel a of Fig. [5] The in- 
elastic relaxation of H 2 (u = l,j = 2) in collisions with 
H 2 (u = 0, j — 0) shown in panel b of Fig. [2] is dominated 
by the transition (1200) — > (1000). In this case, the en- 
ergy gap and the difference in total rotational angular 
momentum of the two CMSs are simultaneously mini- 
mized. The inelastic relaxation of H 2 (v — l,j = 0) in col- 
lisions with H 2 (w = 0,j = 2) is however highly efficient. 
It is dominated by the single transition (1002) — > (1200), 
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FIG. 2: State-to-state cross sections as functions of the final 
combined molecular state for the same systems as in Fig. [1] 
and for a collision energy of 10 -6 K. 



as shown in panel c of Fig. [2] This process conserves the 
total rotational angular momentum and has a very small 
energy gap of 25.45 K. The energy gap arises from the 
slightly different centrifugal distortion of the vibrational 
states v = and v — 1. Therefore, the conservation of 
internal energy and total rotational angular momentum 
is largely satisfied and the process becomes highly effi- 
cient and selective. This near-resonant transition is 150 
times more efficient than the transition (1200) — » (1000). 
The reverse transition (1200) — * (1002) is responsible 
for the rapid enhancement of the inelastic cross section 
for H2(w = l,j = 2) + H 2 (v = Q, j = 0) collisions at 
25.45 K shown in Fig. [1] (the state (1002) becomes ener- 
getically accessible at this collision energy). The mech- 
anism is thus a highly specific process, reminiscent of 
quasi-resonant energy transfer in collisions of rotation- 
ally excited diatomic molecules with atoms [HI HE H(| , 
but with a purely quantum origin. 

In order to elucidate the mechanism for the (1002) 
— > (1200) transition, we have repeated the calculation 
in Fig. [2] for the collision process H 2 (i> = l,j = 0) + 
= 0,j = 2) assuming that the colliding molecules 
are distinguishable. The inelastic transition (1002) — » 
(1200) in collisions of indistinguishable molecules corre- 
sponds to two transitions: [10;02] -> [12;0G] and [10;02] 
— > [00; 12] in collisions of distinguishable molecules. The 
first one corresponds to purely rotational energy transfer, 



and the second to purely vibrational energy transfer. As 
follows from Fig. [2J the rotational energy transfer domi- 
nates over the vibrational energy transfer by three orders 
of magnitude. The inelastic transition (1002) — > (1200) 
can thus be viewed mainly as the rotational energy trans- 
fer. These symmetry considerations can also explain why 
Panda et al. [13] found a surprisingly small cross section 
for the resonant distinguishable process p-H2(v = 0,j = 
0) + o-U 2 {v = l,j = 1) -> p-H 2 (w = l,j = 0) + o- 
H 2 (w = 0, j = 1), despite a small energy gap of 8.5 K. 
The rotational energy transfer between j = and j ' = 1 
is forbidden due to symmetry so this transition is driven 
by the vibrational energy transfer. We conclude that the 
near-resonant transitions should be generally smaller in 
ortho-para collisions than in para-para collisions. 
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FIG. 3: Inelastic cross sections as functions of the initial com- 
bined molecular state for a collision energy of 10 -6 K. 



Fig. [3] demonstrates the effects of vibrational excita- 
tion of the colliding molecules on the inelastic relaxation 
cross section at a collision energy of 10~ 6 K. The results 
are presented as functions of the initial CMSs (uOt/0) 
(dashed line), (u2u'0) (dashed bold line) and (vOv'2) 
(solid line). We use the full basis set for the calcula- 
tion involving the initial CMSs (QjOf), (ljOf), (2j0j') 
and (ljl£) and truncated the basis set, as described 
in Ref. [33|, for the calculations focusing on higher vi- 
brational states. The truncation eliminates all channels 
corresponding to Av and/or Av' > 1 transitions. The 
magnitude of the inelastic cross sections for collisions of 
H2(t>, j = 0) with H 2 (w',j = 0) depends strongly on the 
initial vibrational states of the colliding molecules which 
is consistent with the previous analysis of vibrational re- 
laxation in H-H2 0|, He-Ff.2 Q, and Li-Li 2 [l(| collisions. 
For v = [1 — 6] and v' = 0, the cross sections show a 
monotonous increase with the vibrational quantum num- 
ber v as observed for He-H 2 collisions Q at ultralow en- 
ergies. When both molecules are initially vibrationally 
excited as in v = [2 — 4] and v' = 1, the inelastic cross 
sections become significantly larger and less sensitive to 
the initial vibrational levels of the two molecules. For 
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collisions of H2(v, j = 2) with M^iv' , j = 0), the inelastic 
cross sections are quite large and show a slight increase 
with the vibrational quantum numbers. Finally, for col- 
lisions of H2(v,j = 0) with H2(i/, j = 2), the inelastic 
cross sections are very large due to the near-resonant 
energy transfer and do not display any particular depen- 
dence on the vibrational levels. We have found that the 
elastic cross sections generally do not depend strongly on 
the vibrational excitation of the molecules. They rather 
depend on quantum statistics effect, whether or not the 
vibrational and rotational quantum numbers of the col- 
liding molecules are the same [38| . This difference can be 
seen in the elastic cross sections of (0000) and (1000) in 

Fig.ru 

In summary, we have presented the first study of vi- 
brational relaxation in cold and ultracold molecule - 
molecule collisions and found a new mechanism of state 
selective ro-vibrational energy transfer. If the colliding 
molecules are in different levels of vibrational and rota- 
tional excitation, the ro-vibrational energy may be effi- 



ciently exchanged between the molecules, leading to near- 
resonant energy transfer. If the molecules are distinguish- 
able, this process corresponds to rotational energy trans- 
fer from one vibrational state of the colliding molecules 
to another. This inelastic process that dominates over all 
other inelastic transitions is largely insensitive to the ini- 
tial vibrational excitation of the molecules. The process 
may be an important mechanism for collisional energy 
transfer in ultracold molecules formed by photoassocia- 
tion of ultracold atoms and for chemical reactions pro- 
ducing identical molecules. We hope that these results 
will stimulate further studies of vibrational energy trans- 
fer in molecular collisions at ultracold temperatures. 
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